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> 1 INTRODUCTION 



ABSTRACT 

We study the regeneration of stellar bars triggered by a tidal interaction, using nu- 
merical simulations of either purely stellar or stellar -|-gas disc galaxies. We find that 
interactions which are sufficiently strong to regenerate the bar in the purely stellar 
models do not lead to a regeneration in the dissipative models, owing to the induced 
gas inflow in those models. In models in which the bar can be regenerated, we find a 
tight correlation between the strength and the pattern speed of the induced bar. This 
relation can be explained by a significant radial redistribution of angular momentum 
in the disc due to the interaction, similar to the processes and correlations found for 
isolated barred spirals. We furthermore show that the regenerated bars show the same 
dynamical properties as their isolated counterparts. 

Key words: galaxies: evolution - galaxies: interactions - galaxies: structure - galax- 
ies: kinematics and dynamics. 



Barred galaxies amount to more than one third of the cat- 
alogued disc galaxies in optical wavelengths (e.g., de Vau- 
coulcurs 1963) and to about two thirds in the near-infrared 
(Eskridge et al. 2000| l. It has been established by now, that 
the life of a bar can be divided into several episodes: its for- 
mation, evolution, dissolution and maybe its regeneration 
(e.g., see review by Friedli 1999). 

Two mechanisms for the formation of bars in disc galax- 
ies are being widely accepted at present. Bar formation can 
occur, as shown in many numerical simulations, sponta- 
neously by a global instability in cold, rotationally supported 
stellar discs (e.g.. Miller, Prendergast & Quirk 1970; Hohl 
1971; Ostriker & Peebles 1973; Sellwood 1981; Athanassoula 
& Sellwood 1986; etc.). The second mechanism, also con- 
firmed in numerous A''-body simulations, is the formation of 
tidally induced bars, triggered by interactions with neigh- 
bouring galaxies (e.g., Byrd et al. 1986; Noguchi 1987, 1988, 
1996; Gerin, Combes & Athanassoula 1990; Barnes & Hern- 
quist 1991; Salo 1991; Miwa & Noguchi 1998; etc.). To dis- 
tinguish between the two bar forming mechanisms, Miwa & 
Noguchi 1199811 compared the properties of spontaneously 
formed and tidally induced bars by means of A^-body sim- 
ulations and argued that tidally induced bars are slow ro- 
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tators, while the spontaneously formed bars have usually 
higher pattern speeds. 

Athanassoula 1200311 argued that the evolution of bars 
in isolated disc galaxies is driven by the redistribution of 
angular momentum. Since this is crucial for the work pre- 
sented here, we will briefly summarise the main results. 
Disc galaxies strive to transfer angular momentum out- 
wards jLynden-Bell fc Kainajs 1972^ . Disc material in the 
inner disc, being at resonance with the bar, emits an- 
gular momentum to resonant material in the outer disc, 
or in the spheroidal components like the halo and the 
bulge lAthanassoula 200211 . The corotation radius of the bar 
roughly divides the regions of disc emitters from that of 
disc absorbers. The size and the pattern speed of the bar 
responds to the angular momentum exchange in a way as 
to keep an equilibrium between the emitters and absorbers. 
These recent results argue strongly that the strength and 
the pattern speed of the bar are determined by the amount 
of angular momentum exchanged. The bar's slowdown rate 
found in numerical simulations depends on the relative halo 
mass and on the velocity dispersion of both the disc and the 
halo. Furthermore, Athanassoula 12UU2t has shown that bar 
growth in the disc can be stimulated even in massive halos 
due to the destabilising influence of resonant stars in the 
halo. 

Another factor influencing the evolution of bars is the 
interstellar medium in the disc. Owing to the gravitational 
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torques of the bar, gas can be driven towards the centre 
of the galaxy and this is hkely to be accompanied by cen- 
tral starbursts and formation/fueling of an active galactic 
nucleus (Shlosman, Frank & Begelman 1989). Fully self- 
consistent numerical simulations including gas have shown 
that substantial gas inflow can weaken significantly, or even 
destroy, the bar (Friedli & Benz 1993; Berentzen et al. 1998). 
The central mass concentrations and the times-scales found 
in these simulations, however, might be too high and too 
short, respectively, to be in agreement with the relatively 
high fraction of barred galaxies observed. Furthermore, this 
fraction seems to be independent of galaxy morphology, i.e. 
the same in early- and late-type disc galaxies, while early- 
types are known to have considerably less gas. 

A regeneration of the stellar bar, i.e. a secondary 
episode of bar-formation during the lifetime of the disc, has 
been suggested as a possible mechanism to explain the ob- 
served number of barred galaxies. One scenario for this (Sell- 
wood & Moore 1999; Bournaud & Combes 2002) could be 
the accretion of gas-rich companion galaxies or freshly in- 
falling gas, which, by adding colder material to the disc, 
may cool it sufficiently and allow the generation of a new 
bar. Another scenario suggested, and upon which we will fo- 
cus in this paper, could be the regeneration of a previously 
dissolved or weakened bar by a tidal encounter with a neigh- 
bouring galaxy. In this work we study this regeneration pro- 
cess by means of fully self-consistent, 3D numerical simula- 
tions. The host galaxy, in which the first, i.e. spontaneously 
formed, bar has significantly weakened before the interac- 
tion, is tidally perturbed by a companion galaxy. Various 
orbits and masses have been considered for the companion, 
in order to cover a wide parameter space. We provide a com- 
parison between pure stellar and two-component stars-|-gas 
models and describe the regeneration process and the dy- 
namical properties of the tidally induced bars. 

In <J21we describe the numerical methods and the initial 
conditions of the galaxy models. The evolution and the dy- 
namical properties of the isolated models with and without 
gas are described in 33] and f01 respectively. In JS] and 
we present the results of the interacting models with and 
without gas, respectively. The results are then discussed in 
f[7|and finally we give a summary in JS] 

2 NUMERICAL METHODS AND INITIAL 
CONDITIONS 

2.1 Methods 

To evolve both the coUisionless component, representing the 
stars and the dark matter, and the dissipative component, 
representing the gas, we use an Ai'-body algorithm, com- 
bined with a smoothed particle hydrodynamics (hereafter 
SPH) algorithm (e.g., review by Monaghan 1992). For this 
purpose we use the version of the hybrid 7V-body/SPH code 
of HeUer (1991; see also Heller & Shlosman 1994). The fully 
self-consistent 3D algorithm employs such features as a spa- 
tially varying smoothing length, a hierarchy of time-bins to 
approximate individual time-steps, a viscosity "switch" to 
reduce the effects of viscous shear, and the special-purpose 
hardware GRAPE-3AF to compute the gravitational forces 
and neighbour interaction lists (Sugimoto et al. 1990; Stein- 
metz 1996). 



Table 1. Initial model parameters. 
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2.2 Initial conditions 

Our isolated model 10 is initially composed of a stellar and a 
gaseous disc, embedded in a spherical halo. The radial sur- 
face density of both the discs follows a truncated Kuzmin- 
Toomre (hereafter KT) profile (Kuzmin 1956; Toomre 1963) 
with a radial scalelength ad. Their vertical density pro- 
file follows the sech'^ distribution of an isothermal sheet 
l |Spitzer 1942} . The halo is set up initially as a Plummer 
(hereafter PI) sphere llPlummer 191111 with a radial scale- 
length bh, and is then allowed to relax in the gradually in- 
troduced potential of the frozen disc. For a detailed descrip- 
tion of the initial model 10 see also Berentzen et al. 1200311 . 
The initial model parameters are summarised in Tab.Q The 
columns from left to right give the component, the type 
of density profile, the particle number, the total mass, the 
radial scalelengths ad and &h, the cut-off radius rcut, and 
the vertical scaleheight zq of each component. Our isolated 
model II differs from model 10 in that about one third of 
the gas particles have been replaced with stellar particles, 
i.e. with keeping their positions and assigning new masses, 
keeping the total disc mass the same as in model 10. The 
remaining gas particles have been completely removed from 
the disc. We then assign velocities to the new stellar parti- 
cles according to the mean velocity distribution of the old 
underlying stellar disc. Both host galaxies, 10 and II, are 
evolved in isolation first. The disc of 10 is constructed so as 
to be globally unstable to non-axisymmetric perturbations 
and form a large-scale bar. Owing to induced gas-inflow, the 
stellar bar in model 10 weakens significantly during its evo- 
lution, but is present in both models before the companion 
galaxy is introduced. 

The companion galaxy in the interaction simulations 
is represented by a smoothed point mass (hereafter pt) for 
simplicity, applying a Plummer softening for the force calcu- 
lation. Since the GRAPE-3 hardware can only handle a fixed 
softening length, we apply the same softening for all parti- 
cles, including the companion. The mass ratio Mhost : Mcomp, 
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Figure 1. Evolution of tlie isolated models. We show the face-on isodensity contours of the stellar disc for models 10 (upper panels) and 
II (lower panels) and a grey-scale density plot of the gas particle distribution for model 10. The disc rotation is counter-clockwise, and 
the disc has been rotated in each frame so that the bar is oriented along the x-axis. We give the time on the top of each panel. Time 
and length scales are both in computer units. 



where Mhost and Mcomp are the total mass of the host and of 
the companion galaxy, respectively, has been varied between 
1 : 1 to 1 : 4. In addition a few simulations have been run with 
less massive companions. For the simulations presented in 
this paper we confine ourselves to planar fly-by encounters, 
i.e. the orbits of the host and the companion galaxy are un- 
bound and Me in the equatorial plane of the host. The initial 
conditions have been chosen such as to provide a strong per- 
turbation to the disc. The orbit of the companion galaxy is 
prograde with respect to the rotation in the host disc and 
has been chosen to be either parabolic or hyperbolic, with 
different pericentric separations Rpsri- The pericentric posi- 
tion has been chosen such that at time of pericentre (here- 
after tperi) the major axis of the relatively weak bar, which 
has formed/evolved in the hosts disc before the encounter, 
points towards the companion. The times of pericentric sep- 
aration have been chosen to be fpori = 270 and 500 for all 
simulations with host galaxy 10 and II, respectively. The 
initial orbital parameters at tpai have been calculated from 
the solution of the corresponding two-body problem, either 
for a given orbital eccentricity e, or for a pericentric angular 
frequency of the companion 

^pcri — ^pcri/ Rpcri , (1) 

where fperi is the velocity of the companion at pericentre. 
To obtain the initial centre of mass positions and velocities 
of the host and of the companion galaxy, we integrate their 
orbits backward in time, starting from the pericentric time 
tperi, until the distance AR between host and companion is 
roughly 10 times the cut-off radius of the initial halo. During 
this backwards orbit integration the particles of the host are 
frozen with respect to each other. Owing to the limited spa- 
tial range of the GRAPE-3 hardware, the interaction models 
have been calculated in the inertial frame of the host galaxy 
and the force of the companion has partly been calculated 
by direct summation on the front-end, i.e. without using the 
GRAPE hardware. 



2.3 Units 

The adopted units for mass, distance, and time are M = 
6 X 10^° M0, R = 3kpc and r = lO'^ yr, respectively, for 
which the gravitational constant G is unity. The dynamical 
time is Tdyn = {rf^/G Mi^)^^'^ = 4.8 x lO'' yr, where Mi^ 
is half the total mass of the host and is its half-mass 
radius, which is, after relaxation of the halo, approximately 
8.5 kpc or in terms of the disc scalelength, 2.8 ad- The initial 
stellar disc rotation period in these units then corresponds 
to trot = 27rrdyn « 3 X 10** yr. A fixed gravitational softening 
length of 6 = 0.375 kpc is used for all particles. An isothermal 
equation of state is used for the gas with a sound speed of 
Ws = 12km s~^. The corresponding thermal gas temperature 
is 10'' K. 



3 ISOLATED MODEL WITH GAS - MODEL 10 

As described in the previous section, the isolated galaxy 
model 10 is composed of a stellar (coUisionless) and a gaseous 
(dissipative) disc component, embedded in a Uve halo. The 
morphological evolution of the disc is shown in Fig. (upper 
panels). The stellar disc has been chosen to be initially bar 
unstable and forms a large-scale bar within roughly 2 trot- 
The stellar bar reaches its maximum strength at about t = 60 
(see Fig. |2l . The bar strength, here defined by the normal- 
ized m = 2 Fourier component of the stellar disc mass dis- 
tribution, has been measured inside a cylindrical radius of 
7?= 1.25 and about one scaleheight of the disc plane. Owing 
to the torque of the bar, most of the gas in the bar region 
is driven towards the centre of the galaxy, accumulating in 
a nuclear and circumnuclear disc. At the end of the run the 
mass of the gaseous nuclear disc represents some 18 per cent 
of the total mass within a spherical radius of 1.0 kpc. As a 
result of the growing central mass concentration the stel- 
lar bar weakens significantly, and the disc settles down in a 
quasi-stable state at about t ~ 120. Some gas also accumu- 
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lates near the end of the weak stellar bar, forming an oval 
ring close to the position of the inner ultra-harmonic reso- 
nance (hereafter UHR) at R — 2.25 kpc. The pattern speed of 
the bar (hereafter Qp) first increases due to the short burst 
of gas inflow and thereafter reaches a constant rate of about 
f2p = 0.3T~^, or 29.3 km s~^ kpc~^. The evolution of model 
10 is followed in total up to t = 300, i.e. about 10 disc ro- 
tations. A more detailed description of the evolution of this 
model and its dynamical properties is given in Berentzen et 
al. (20031. 

To identify the presence and the location of the main 
planar resonances in the disc, we apply both the lin- 
ear (epicyclic) approximation and non-linear methods. We 
therefore calculate the gravitational potential in the z = 
plane of the disc on a Cartesian grid and symmetrise it with 
a four-fold symmetry with respect to the main axes of the 
bar. From the potential we then derive the azimuthally av- 
eraged circular and epicyclic frequency Qc{R) and n{R), re- 
spectively. The standard linear resonance condition is given 
by 



Qc -\ K 

m 



(2) 



with integers I and m. For a strict definition of the in- 
ner Lindblad resonances (hereafter ILRs; with I = —1 and 
m = 2) three extensions to the linear definition have been 
proposed (see Athanassoula 1200311 for a discussion), which 
are not fully compatible to each other and sometimes in con- 
tradiction. Throughout this work, unless otherwise noted, 
we apply the orbital structure definition for the ILRs, i.e. 
by saying that an ILR exists if and only if the X2 and 
orbit families exist. The ILRs have been confirmed in this 
case by calculating the surfaces of section (hereafter SOS). 
The SOS have been constructed by integrating orbits of a 
given Jacobian energy (hereafter Ej; see also Sec.|lJ in the 
equatorial plane of the disc and marking the points in the 
(y, y) plane each time the orbits cross the line x = with 
i < (e.g., Binney & Tremaine 1987). At the end of the 
run {t = 300) the corotation radius (hereafter CR; i.e., the 
radius at which f2c = np) in this model is at about 6.15 kpc 
{Ej « —0.91), and the ILRs are located at about 0.47 kpc (in- 
ner ILR; Ej = -1.58) and 2.1 kpc (outer ILR; _Bj = -1.10), 
respectively. 



4 ISOLATED MODEL WITHOUT GAS - 
MODEL II 

We study the regeneration of the bar by interactions of two 
different types of host galaxies - with and without gas in 
the disc - which, nevertheless, have similar dynamical prop- 
erties. We therefore construct our purely stellar model II by 
replacing the gas particles in model 10 with stellar particles, 
as described in Sec. 12.21 This has been accomplished for the 
last snapshot of model 10 at t = 300, i.e. at the end of the 
run, when the amplitude of the initial bar has very consider- 
ably decreased and the disc is in a quasi-stable state. After 
replacing the gaseous particles with stellar ones we follow 
the evolution of the new model II for another At — 300. The 
morphological evolution of the stellar disc is shown in Fig.Q 
(lower panels) , and in Fig. |5| we show the evolution of the 
bar strength and of the pattern speed. The bar strength first 
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Figure 2. Evolution of the bar strength (upper panel) and pat- 
tern speed (lower panel) for models 10 (stars-|-gas) and II (stars). 
The transition between the two models is marked by two vertical 
dashes (t = 300). 




Figure 3. Characteristic diagram of the isolated model II at 
the end of the run. We show the characteristic curves of the main 
orbit families (full lines) and the zero velocity curve (ZVC; dashed 
line). 



decreases abruptly, owing to the newly added stellar parti- 
cles, which follow initially the distribution of the gas. The 
phase angle of the m = 2 Fourier component of the latter is 
shifted with respect to that of the old underlying disc and 
therefore tends to decrease the normalized m = 2 amplitude 
of the combined discs. The velocities, which are assigned to 
the new stellar disc component, follow the velocity distribu- 
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Figure 4. Examples of the morphological evolution of interaction models with and without gas. We show the face-on isodensity contours 
of the stellar disc for the models I0C4p (with gas, upper panels) and II C4p (without gas, lower panels). For the first model a grey-scale 
plot of the gas particle distribution is shown. The layout is as for Fig. 
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Figure 5. Evolution of the bar strength as a function of time 
for prograde, parabolic encounters with host galaxy 10 (thick 
line). The pericentric separation is ijp(,i-i = 15 and the mass ratio 
-'^^host ■ -'\<^comp is given on top of each frame. The type of the com- 
panions orbit is given in the upper right corner of each frame. For 
comparison we also show the bar strength of the isolated model 
(thin line). 



tion of the old underlying stellar disc. Therefore the stellar 
disc of model II is initially slightly out of virial equilibrium, 
since the dispersion of the stellar particles is higher than that 
of the gas particles. After some At = 30, however, the disc 
settles to a new equilibrium and the bar strength increases 
again, reaching its previous value and remaining constant 
after the short relaxation phase. Owing to the decreasing 
central mass concentration of the new stellar component, 
the pattern speed of the bar slightly decreases first, but re- 
mains roughly constant after the relaxation phase at about 
f2p = 0.22r~^, or 21.5 km s~^ kpc~^, till the end of the run. 

To gain some insight in the orbital structure of the disc 
of model II, we locate the main planar periodic orbits in 
a frozen potential that rotates with the bar. General infor- 



mass ratio 1 :4 

R.„-:=15 



pen 



^ 0.3 



P2 - 




t [lO^yr] 



t [lO^yr] 



t [lO^yr] 



Figure 6. Same as Fig. |^ but for encounters with a fixed mass 
ratio of 1 : 4 and pericentric separations of 15 (p), 10 (p2) and 5 
(pa), respectively. 



mation on orbits in barred potentials can be found, e.g., in 
Contopoulos & Grosb0l (1989). The gravitational potential 
has been calculated the same way as described in the pre- 
vious section, but is time-averaged over roughly one bar ro- 
tation. For simplicity, we restrict ourselves to planar orbits 
which are bi-symmetric with respect to the bar and close 
after one orbit around the centre in the frame of reference 
corotating with the bar. The results are then displayed in 
terms of a characteristic diagram, where for each orbit we 
plot its Jacobi integral E} with respect to the j/-intercept 
value with the j: = plane. The Jacobi integral is a con- 
served quantity along any given orbit in the rotating frame, 
and can be thought of as an effective energy (e.g., Binney 
& Tremaine 1987). In the characteristic diagram the orbits 
form curves of families. In Fig. I^we show the characteris- 
tic diagram of the main orbit families, using the notation 
of Contopoulos & Papayannopoulos (19801, in model II at 
the end of the run (t = 600). The dashed curve is the zero- 
velocity curve (hereafter ZVC), which delineates the acces- 
sible region in the plane based on energy considerations. 
The family labelled xi consists of orbits that are elongated 
along the bar and predominately gives the bar its structure. 
Orbits of the X2-family are elongated perpendicular to the 
bar. Their presence is indicative of an ILR (or more than 
one) in the non-linear regime. The X4 orbits are retrograde 
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and slightly elongated perpendicular to the bar. The coro- 
tation radius in model II is at Ej = —0.79, or approximately 
7.98 kpc. The j:2-orbits range in energy from = — 1.35 to 
— 1.07, with semi-major axes ranging from 0.86 to 2.1 kpc. 



5 INTERACTING MODELS WITH GAS 

In this section we describe the simulations in which the gas- 
rich host galaxy 10 is perturbed by encounters with compan- 
ions of different mass. The time of pericentric separation has 
been chosen to be tperi = 270 for all runs with the host galaxy 
10. In the first set of simulations the companions are initially 
set up to follow prograde, parabolic orbits with a pericen- 
tric separation of Rpai = 15, which is just outside the halo of 
the host galaxy. When referring to a specific model, we will 
hereafter use a notation like, for instance, 'I0C4p', where 
'10' is the name of the host galaxy, 'C4' is the name of the 
companion (see Tab.0, and 'p' (parabolic) denotes the type 
of orbit, the latter following the notation given in Tab.|21 In 
Tab.|21we summarise the basic parameters and properties of 
both the isolated and the interaction models. The first and 
the second column give the name of the host and of the com- 
panion galaxy, respectively, and the third column gives the 
orbit of the companion (see Tab.j^J. In the fourth column we 
give the time when the companion is included in the simu- 
lation, resulting from the initial backward orbit integration 
(see Sec. 12.21 . The fifth to seventh column give the peri- 
centric separation and frequency, and the eccentricity of the 
orbit, respectively. In the eighth and ninth column we give 
the strength and pattern speed of the bar, respectively, after 
the interaction. The tenth and eleventh columns give the ra- 
tio of bar strength and pattern speed to the corresponding 
values of the isolated model. In the penultimate column we 
give the interaction strength parameter as defined in Sec.|S]l 
and in the last column the corotation radius at the end of 
each run, determined from the linear analysis. 

In Fig. 0] (upper panels) we show an example of the 
morphological evolution of the disc during an encounter, in 
this case for model I0C4p. The bar strength as a function 
of time for the models with different massive companions 
is shown in Fig. |K| We find that the encounters in this set 
of simulations are not sufficiently strong to regenerate the 
stellar bar in the disc, even for a mass ratio of 1 :4, i.e. even 
with companion C4. In order to increase the strength of the 
tidal perturbation on the host galaxy, we run two additional 
simulations with companion C4, in which we now choose 
pericentric separations of 7?pcri ~ 10 and 5. The orbits of the 
companion are again chosen to be prograde and parabolic. 
The bar strength for these models is shown in Fig. |n] The 
peaks immediately following tpcri result from the time when 
the transient tidal arms contribute to the m — 2 power in the 
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Figure 7. Logarithmic surface density as a function of radius 
for the stellar and gaseous disc and of the halo (from top to bot- 
tom, respectively). We show logS(/?) for models 10 C4p (left-hand 
panels) and II C4p (right-hand panels). The dotted and full lines 
represent times before (At = —100) and after (At = 240) the in- 
teraction, respectively. 
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Figure 8. Specific cumulative angular momentum as a func- 
tion of radius in the stellar disc at the end of the run of the 
isolated model ID (with gas) and of the interaction models with 
a mass ratio of 1 :4 and different pericentric separations. The full 
and dashed lines show the results for the stellar and gaseous disc, 
respectively. 

inner disc region, in which the bar amplitude is measured, 
but disappear as soon as the tidal arm features have dis- 
solved. Even with these close encounters, it is not possible 
to regenerate the stellar bar in the disc of the gas-rich host 
galaxy. 

In order to facilitate comparison with the purely stellar 
simulations described in the next section, we present here 
a more detailed analysis of the interaction simulations with 
the gas-rich host galaxy 10. As can be seen in Fig. |1] already, 
the interaction apparently leads to a significant redistribu- 
tion of both the stellar and the gaseous material in the disc. 
This becomes more evident in Fig. |7| in which we show the 
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Figure 9. Same as Fig.|^ but for the combined disc (stars+gas). 
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Figure 10. Bar strength as a function of time for the prograde, 
paraboUc encounters with the purely stellar model II (thick line). 
The layout is as in Fig. ISl 



logarithmic radial surface density of the disc and the halo 
for model I0C4p at times before and after the interaction. 
While the stellar density remains basically constant within 
a radius of approximately J?<2.0, i.e. roughly inside corota- 
tion in the corresponding isolated model, there is some net 
inflow of gas from the outer disc. The inflowing gas first piles 
up at a radius of about 1.25, which is close to the po- 
sition of the gaseous ring, and then flows gradually towards 
the central disc region, accumulating in the nuclear disc. At 
the end of the run, the gas mass within 0.5 length units has 
increased by about 0.5 per cent of the total host mass, or 
some 6 per cent of the total gas mass, while the gas mass 
within 2 length units increased by roughly 1 per cent of the 
total mass, or about 12 per cent of the total gas mass. Some 
fraction of both the stellar and the gaseous material in the 
outer disc moves further out, contributing to the transient 
tidal arms and/or getting stripped off by the companion. 
The halo material in the inner region does not show any 
net radial redistribution, but we find some expansion in the 
outer parts, where halo material is also stripped off by the 
companion. 

We measure also the specific cumulative angular mo- 
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Figure 11. Characteristic diagram of model IlC4p after the 
interaction. The layout is as in Fig. |3] 



mentum (hereafter L^) as a function of radius for both the 
stellar and the gaseous disc, given by 




(3) 



where Ri, rrn and Zz,i are the planar radius, the mass and the 
angular momentum, respectively, of the i-th disc particle. 
As an example we show in Fig. |H|the results of simulations 
with companion C4 and different pericentric separations. As 
can be seen in this plot, the interaction removes angular 
momentum from the stellar disc at all radii shown. We do not 
find any indication that the corotation radius of the bar or 
of the companion separates regions which gain from regions 
which lose angular momentum. More angular momentum 
is removed from the stellar disc for the closer pericentric 
separation. The gaseous disc shows a different behaviour. 
The gas in the inner disc gains angular momentum compared 
to the isolated model, within roughly R = 3.0 and 1.7 for 
pericentric separations of 7?peri = 15 and 10, respectively, 
and loses angular momentum outside those radii. 

In Fig. |5] we show the specific cumulative angular 
momentum as a function of radius in the combined disc 
(stars+gas) . We find that the disc in the interaction models 
10 C4 p and 10 C4 p2 does not lose, or even gains, angular mo- 
mentum within a radius of _R~2.0 and i?«1.0, respectively, 
as compared to the isolated model 10. 



6 INTERACTING MODELS WITHOUT GAS 

In this section we describe the set of simulations with 
the purely stellar model II as the host galaxy. We start 
with simulations having similar orbits for the interaction 
as the ones described in the previous section, i.e. prograde 
and parabolic encounters using different companions (see 
Tab. Ol. As shown in Fig. 1101 in which we plot the bar 
strength as a function of time, simulations with mass ra- 
tios of 1 :3 and 1 :4 are, contrary to the dissipative models, 
sufficiently strong to regenerate the stellar bar in the disc. 
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Table 3. Main model parameters and properties. 



Host 


Comp 


orbit 


^start 


^pcri 


'^pcri 


e 




Qp 








CR 


10 


" 




0.0 








0.14 


0.30 


1.0 


1.0 




2.05 


11 






100.0 








0.16 


0.22 


1.0 


1.0 




2.76 


10 


CI 


p 


170.0 


15.0 


0.049 


1.00 


0.10 


0.30 


0.73 


1.33 


1.927 


1.80 


10 


C2 


p 


170.0 


15.0 


0.060 


1.00 


0.13 


0.30 


0.91 


1.01 


3.159 


2.13 


10 


C3 


p 


170.0 


15.0 


0.069 


1.00 


0.08 


0.27 


0.61 


0.92 


4.116 


2.31 


10 


C4 


p 


170.0 


15.0 


0.077 


1.00 


0.11 


0.26 


0.77 


0.87 


4.922 


2.47 






P2 


170.0 


10.0 


0.141 


1.00 


0.16 


0.18 


1.19 


0.63 




2.99 






P3 


170.0 


5.0 


0.400 


1.00 


0.10 


0.22 


0.70 


0.74 






10 


C4 




250.0 


10.0 


0.300 


8.00 


0.06 


0.28 


0.44 


0.93 


1.906 


2.20 


11 


C5 


g 


450.0 


5.0 


0.300 


3.22 


0.18 


0.20 


1.14 


0.90 


0.717 


2.94 


11 


C6 


J 


470.0 


7.5 


0.300 


11.66 


0.17 


0.21 


1.08 


0.92 


0.419 


2.94 


11 


C6 


s 


450.0 


5.0 


300 


2 75 


0.21 


0.19 


1 31 


0.82 


1.105 


2 78 


11 


CI 


p 


400.0 


15.0 


0.049 


1.00 


0.16 


0.22 


1.00 


0.98 


1.927 


2.87 


11 


CI 


(J 


480.0 


10 


300 


21.50 


0.17 


0.21 


1 11 


0.93 


0.458 


2 95 


11 


CI 




490.0 


10.0 


0.600 


89.00 


0.16 


0.21 


1.04 


0.95 


0.261 


2.90 


11 


CI 


g 


440.0 


5.0 


0.300 


1.81 


0.27 


0.12 


1.74 


0.56 


2.417 


4.13 


11 


CI 




480.0 


5.0 


0.600 


10.25 


0.22 


0.18 


1.38 


0.82 


1.124 


3.25 


11 


CI 




480.0 


5.0 


0.900 


24.31 


0.19 


0.19 


1 19 


0.87 


0.833 


3 20 


11 


CI 


J 


470.0 


7.5 


0.300 


8.48 


0.21 


0.19 


1.33 


0.84 


0.852 


3.19 


11 


C2 


p 


400.0 


15.0 


0.060 


1.00 


0.17 


0.21 


1.08 


0.94 


3.159 


2.95 


11 


C2 


d 


480.0 


10.0 


0.300 


14.00 


0.18 


0.20 


1.15 


0.92 


0.928 


2.99 


11 


C2 




490.0 


10.0 


0.600 


59.00 


0.17 


0.21 


1.09 


0.95 


0.523 


2.90 


11 


C2 




480.0 


5.0 


600 


6.50 


0.21 


0.18 


1 33 


0.81 


2.303 


3 27 


11 


C2 




480.0 


5.0 


900 


15.88 


0.20 


0.19 


1 25 


0.87 


1.679 


3 10 


11 


C3 


p 


400.0 


15.0 


0.069 


1.00 


0.21 


0.18 


1.36 


0.81 


4.116 


3.30 


11 


C3 


a 


480.0 


15.0 


0.300 


36.97 


0.17 


0.21 


1.05 


0.94 


0.500 


2.86 


11 


C3 


d 


480.0 


10.0 


0.300 


10.25 


0.21 


0.18 


1.36 


0.82 


1.410 


3.29 


11 


C3 


e 


490.0 


10.0 


0.600 


44.00 


0.19 


0.20 


1.18 


0.89 


0.787 


3.07 


11 


C3 


f 


490.0 


10.0 


0.900 


100.25 


0.17 


0.20 


1.09 


0.94 


0.594 


2.91 


11 


C3 


h 


480.0 


5.0 


0.600 


4.63 


0.26 


0.13 


1.64 


0.61 


3.545 


3.71 


11 


C3 


i 


480.0 


5.0 


0.900 


11.66 


0.24 


0.17 


1.51 


0.74 


2.541 


3.35 


11 


C4 


P 


400.0 


15.0 


0.077 


1.00 


0.23 


0.17 


1.45 


0.77 


4.922 


3.43 


11 


C4 


a 


480.0 


15.0 


0.300 


29.38 


0.17 


0.21 


1.09 


0.94 


0.803 


2.93 


11 


C4 


b 


490.0 


15.0 


0.600 


120.50 


0.16 


0.21 


1.05 


0.95 


0.460 


2.88 


11 


C4 


c 


490.0 


15.0 


0.900 


272.38 


0.16 


0.22 


1.03 


0.97 


0.349 


2.82 


11 


C4 


d 


480.0 


10.0 


0.300 


8.00 


0.24 


0.17 


1.52 


0.74 


1.906 


3.51 


11 


C4 


e 


490.0 


10.0 


0.600 


35.00 


0.17 


0.20 


1.11 


0.90 


1.052 


3.05 


11 


C4 


f 


490.0 


10.0 


0.900 


80.00 


0.17 


0.21 


1.10 


0.93 


0.793 


2.94 


11 


C4 


h 


470.0 


5.0 


0.600 


3.50 


0.22 


0.13 


1.41 


0.57 


4.869 


3.66 


11 


C4 


i 


480.0 


5.0 


0.900 


9.13 


0.24 


0.16 


1.50 


0.70 


3.419 


3.41 



and the strength of the induced bar increases with the mass 
of the companion. In Fig. 1111 we show as an example the 
characteristic diagram of the host galaxy of the interaction 
model II C4 p at time t = 640, i.e. the end of the run, in which 
the stellar bar has been successfully regenerated by the in- 
teraction. The orbit analysis has been performed following 
the description given in Sec.|l] The corotation radius in this 
model is at Ej = —0.63, or about 10.3 kpc. The X2 orbits 
range in energy from = — 1.34 to —0.88, with semi-major 
axes ranging from 0.83 to 3.04 kpc. 

Motivated by these results, we now explore a much 
wider - though still far from complete - parameter range, in 
order to determine the initial conditions necessary (or suffi- 
cient) to regenerate the stellar bar in the disc. We therefore 
run a set of simulations, in which we vary the orbit of the 



companion, the pericentric separation 7?pori and the mass 
Mc , respectively. The simulations and their basic initial pa- 
rameters are summarised in Tab. |3| From the solution of 
the corresponding two-body problem the eccentricity e of 
the companions orbit has been determined by its angular 
frequency Upai at pericentre, and we restrict ourselves in 
this work to unbound orbits only. We choose Li^peri such that 
a prescribed number of inner Lindblad resonances is present 
in the disc, i.e. that certain orbits in the disc are in reso- 
nance with the companion, when passing Rpcri-^ Since the 
stellar bar in the unperturbed host galaxy is very weak, we 

^ We point out that the angular frequency of the companion, and 
therefore also the location of the corresponding resonances in the 
disc, changes with time. 
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Figure 12. Basic frequencies of tlie isolated model II (t = 500) 
as a function of radius, obtained with the hnear axisymmetric 
definition, using the epicyclic approximation. The horizontal lines 
mark the angular frequency cupai of the companion at the time 
of pericentric separation. 
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Figure 13. Different trajectories of the companion galaxy C4. 
We show the parabolic and hyperbolic orbits (full lines) from 
simulations with difl^erent pericentric separations, as given on top 
of each panel. The corresponding circular and parabolic orbits are 
indicated by dotted lines. The host galaxy is marked by a star in 
each panel. 



apply the linear approximation, as described in Sec^l in or- 
der to determine the resonances in the disc. The resonance 
diagram for the isolated model II is shown in Fig. 1121 and 
the frequencies chosen for the companion are: 



OJperi — 0.3 
<.^pGri = 0.9 



(two ILRs) 
(one ILR) 
(no ILR) . 



The mass of the companion for these simulations has been 
varied, as before, from 1 : 1 to 1:4, and we have a few ad- 
ditional runs with even lower mass ratios, namely 1 : i (C6) 
and 1 : i (C5). A complete list of all simulations and of their 
main parameters is given in Tab.|3 Some of the companion 
orbits are illustrated in Fig. 1131 for different eccentricities 
and pericentric separations. 

In Fig. 1141 we show the bar strength as a function of 
time for this set of interaction models. It is striking that 



- in contrast to the models with gas - the interaction is 
sufficiently strong to regenerate the bar in about half our 
models. As can also be seen from the plot, there is a general 
trend, that the strength of the induced bar increases with 

• increasing mass Mc of the companion, 

• decreasing pericentric separation -Rperi, 

• decreasing pericentric frequency cjpcri- 

This is what one might generally expect, since the gravita- 
tional force of the companion is proportional to its mass, 
the tidal force decreases with R~^, and the time integrated 
force of the companion on the host galaxy becomes stronger 
for slower passages. 

Furthermore, after the bar has formed in the disc, we 
find that the amplitude of the bar stays constant with time. 
We have checked the life-time of the regenerated bar by run- 
ning model II C4 d, which shows a significant increase in bar 
strength after the interaction, for roughly At = 4 x 10^° yr 
in total. Using a linear-fitting we find a decay-rate of the 
bar of roughly 9 x 10^^^ yr. With this the bar will drop to 
half its amplitude after approximately At^SxlO^^yr. After 
one Hubble time (with Hg^ = 1.3 x 10^" yr) the bar ampli- 
tude would have dropped only by 12 per cent. A fraction of 
this decrease could be introduced by the relatively low num- 
ber of particles used in the simulations, so that in fact there 
could be (almost) no sign of a decay for bars formed in inter- 
actions. These results strongly argue for the fact that bars 
formed by the interaction are long-living and by no means 
transient phenomena. 

6.1 Dynamical properties of the bars 

To quantify the correlations described in the previous sec- 
tion, we define a parameter Q, which allows us to evaluate 
the interaction strength: 



e ; 



gal 



./gal 



dt, 



(4) 



where /gai and /comp are the mean radial forces, exerted from 
the host II and the companion galaxy, respectively, averaged 
over the equatorial plane of the host galaxy. In practice we 
calculate them on a radial equally-spaced polar grid with 
maximum radius R — 4.8 and then take a density weighted 
average. The integral in eq. (0J is carried out from t = — 250 
to t = 250, with t = corresponding to tpai, in order to 
guarantee an adequate convergence of O in all models. For 
these force calculations both galaxies are approximated by 
point masses for simplicity. The results for O are given in 
the penultimate column of Tab. |21 We also calculate the 
normalized bar strength A /An, where A and An are the 
bar strength in the interaction model and in the isolated 
model, respectively. By plotting the normalized bar strength 
versus the interaction strength parameter O (Fig. I15II . we 
find a roughly linear correlation between the two quantities 
for each companion Ci, with i — 1 ... 4:. One of the 1 : 4 
models (IlC4h) is systematically off-set in all plots of this 
kind (i.e.. Figures ITSl 1171 and 1 191 . but always lies within the 
2(7 confidence limit. The models with parabolic encounters 
are all off-set from the relationships found, but seem to be 
correlated themselves linearly. 

Plotting the bar strength of the different models in a 
diagram showing the logarithmic orbital eccentricity e versus 
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Figure 14. Bar strength as a function of time for some of the pure stellar models. The four frames show the models with mass ratios 
of 1:1, 1:2, 1:3 and 1:4, respectively. The filled squares mark models which have not been run, because other runs allow us to deduce 
that a bar cannot form. Models in which the companion would be on a bound elliptical orbit are marked by a filled circle. 



the angular frequency LOpai (see Fig. ll6ll . we find that models 
with orbits of roughly log e > 1.2 are not sufficiently strong 
to regenerate the bar in the purely stellar models. A larger 
sample of simulations has to be run, however, to confirm 
the existence of such a limiting eccentricity. If such a limit 
really exists, it would give a necessary, but not sufficient 
condition for tidally induced bar formation for our purely 
stellar model. 

Angular momentum exchange in isolated galaxies leads 
to correlations between the bar pattern speed and strength, 
as well as between the bar strength and the angu- 
lar momentum gained by the spheroid, which, in many 
cases, is a measure of the angular momentum exchanged 
JAthanassoula 200311 . We will now test whether such cor- 
relations can be found in our regenerated bars. We indeed 
find a tight correlation between the strength and the pat- 
tern speed of the regenerated bars (Fig. 1171 . The stronger 
the bar gets after the interaction, the lower becomes its pat- 
tern speed, in good agreement with Athanassoula 121)1)311 . 
We further find that the pattern speed in the interaction 
models is always lower than in the corresponding isolated 
case. 

The passage of the companion again is accompanied 
by a redistribution of disc angular momentum. The angular 
momentum exchange between the different disc regions (as 
defined in the previous section) is similar to that in Fig. |H| 
In Fig. 1181 we show the radial distribution of the specific cu- 



mulative angular momentum L^{R) for a set of simulations 
with companion C4 in comparison with the isolated model 
II. For the runs with lower mass companions, we find similar 
results. It is noticeable that the regeneration of the stellar 
bar has been feasible only in those models which show a 
significant change in Lz- Actually we find a tight correla- 
tion between the loss of angular momentum ALz in the disc 
measured inside the initial disc cut-off radius r-cut and the 
induced bar strength, as shown in Fig. 1191 The more angu- 
lar momentum is removed from the disc, the stronger the 
regenerated bar becomes, in good agreement with what was 
found for isolated bars by Athanassoula II2003II . We find the 
same correlation about equally strong when we plot the an- 
gular momentum change within the corotation radius. This 
important connection between the angular momentum ex- 
change and the (re)generation process of the bar is described 
further in the discussion section. 



7 DISCUSSION 

The regeneration of stellar bars triggered by galaxy interac- 
tions has been suggested as an additional scenario for the 
formation of bars and for explaining the observed frequency 
of bars along the Hubble sequence (Sellwood & Moore 1999 
; Friedli 1999; Athanassoula 2000). So far only the formation 
of tidally induced bars in initially non-barred disc galaxies 
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Figure 15. Normalized bar strengtli A/An versus tlie interaction 
strengtli parameter ©. We sliow tlie results of the purely stellar 
simulations and the corresponding linear fits for mass ratios of 
1:1 (bullets, full line), 1:2 (squares, dashed line), 1:3 (triangles, 
dotted-dashed line) and 1:4 (stars, dotted line). The simulations 
with parabolic orbits are shown with open symbols and are not 
taken into account for the linear fits. 
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Figure 16. Logarithmic orbital eccentricity versus the angular 
frequency of the companion at pericentre time. For each peri- 
centric separation the curves obtained from the corresponding 
two-body problem are plotted for different mass ratios: 1 : 1 (full 
line) ,1:2 (dashed) ,1:3 (dotted-dashed) and 1 : 4 (dotted line) . 
The filled circles mark the different simulations and their size is 
proportional to the bar strength at the end of each run. The two 
circles which do not lie on a line are simulations with mass ratios 
of 1 : 0.5 and 1 : 0.33. The limit of e, as explained in the text, is 
marked by the horizontal dashed line. 

and the required conditions, under which such an event may 
occur, have widely been studied by means of numerical sim- 
ulations (Byrd et al. 1986; Noguchi 1988; Salo 1991; Miwa 
& Noguchi 1998; etc.). These results, however, need not nec- 
essarily apply to interactions with a formerly barred galaxy. 
Indeed, a former stellar bar likely might have changed the 
dynamics of the disc, i.e. increased the velocity dispersion 
in the stellar disc and changed both the density and angular 
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Figure 17. Normalized pattern speed Qp/Qu versus normalized 
bar strength A/ An of the purely stellar interaction models. Sym- 
bols are as in Fig. 1151 
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Figure 18. Specific cumulative angular momentum of the stellar 
disc for the isolated model II (thick lino) and different interac- 
tion models II C4 (thin lines). The dotted lines indicate models 
in which no or only a weak increase of bar strength has been 
found after the interaction. Models in which the bar has been re- 
generated (full lines) are labelled (with the orbit type as given in 
Tab. 13 and the corresponding bar pattern speed is given. 

momentum distribution (compare Fig. 1201 . Therefore, the 
regeneration of stellar bars is subject to different conditions 
than the formation of a bar in an bar-unstable isolated disc 
or in a tidal interaction. 

In this work we present numerical simulations starting 
with a host galaxy which has initially been bar unstable and 
in which the bar has been significantly weakened due to gas 
inflow before the interaction. 

7.1 Quantifying the interaction strength 

To quantify the strength of the interaction, we have intro- 
duced a parameter Q for the interaction strength (Sec. 6.1), 
which basically takes the following quantities into account: 
1) the mass of the companion, 2) the pericentric separation 
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Figure 19. Total change of angular momentum ALz of the stellar 
disc inside its initial cut-off radius as a function of bar strength 
normalized by the corresponding value of the isolated model. The 
different symbols represent interactions with different mass ratios. 
Parabolic and hyperbolic orbits are represented by open and filled 
symbols, respectively. 
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Figure 20. Logarithmic surface density of the stellar disc of mod- 
els 10 (left-hand panel) and II (right-hand panel). The filled dots 
show the density distribution of the stellar discs and the thick 
full lines the best-fit of the corresponding Kuzmin-Toomre pro- 
file. The thin full line in the left-hand panel gives the density 
distribution of the initial stellar disc in model 10. We also show 
the different exponential discs used by Miwa & Noguchi II1998I 
with dotted lines and give the mass of each disc in their model 
units. The value of the Toomre Q-parameter of the stellar disc is 
given at the bottom of each panel. 



between the galaxies, and 3) the velocity of the compan- 
ion at pericentre, defining a kind of interaction time-scale. 
With our definition (see Sec. 6.1), we find that the interac- 
tion strength correlates well, i.e. roughly linearly, with the 
strength of the regenerated bars, which have been created 
in our purely stellar models. 

A different parameter for the interaction strength has 
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Figure 21. Normalized bar strength A/An versus the Dahari in- 
dex Qd (Dahari 1984) for interactions with different mass ratios. 
The layout is as in Fig. 1151 



been introduced by Dahari 119841 as an estimate of the di- 
rect tidal impulse: 



Qd = (Mcomp/Mhost)/(i?miii/i?disc) 



(5) 



where Mc, 



and i?pori are the mass of the companion 



galaxy and the pericentric separation, respectively. Mhoat 
is the mass of the host galaxy within the disc truncation 
radius -Rdisc. Salo (,1,991 ) has used the Dahari index as a 
quantitative measure for tidally induced bar formation in 
his 2D simulations and found a minimum value of Qd, de- 
pending on the specific host galaxy, above which a bar is 
formed by the interaction. The specific limiting value of Qd 
Salo found depends strongly the central mass concentration 
of the host galaxy and cannot be easily transferred to our 
model. Furthermore, applying the Dahari index to our sim- 
ulations we find neither a correlation between the strength 
of the interaction and the regeneration of the bar, nor a 
limiting value of Qd (see Fig. I2H . The main drawback of 
the Dahari index is that it does not take into account the 
interaction orbit, i.e. the interaction time-scale. The Dahari 
index therefore does not seem to be sufficient to constrain 
the parameters necessary for the formation or regeneration 
of bars by galaxy interactions in general. 

A more advanced parameter has been introduced by 
Elmegreen et al. 1199111 . This is based on the Dahari index, 
but includes also the ratio between an interaction time-scale 
AT and some dynamical time T of the disc: 



= Qd X 



AT 
— , 



(6) 



where T = i?gai/GA/gai and AT is the time it takes the com- 
panion to move by one radian relative to the hosts centre at 
pericentre time. These authors also report a limiting value 
Qe = 0.038 for the formation of the bar by tidal interactions. 
Applying this parameter to our simulations we find a cor- 
relation between the interaction strength and the strength 
of the regenerated bar (see Fig. I22II . This, however, is less 
pronounced than the corresponding one for our parameter 
Q, which takes more fully the interaction into account. The 
existence of a limiting parameter, as found by Salo 119911 
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Figure 22. Normalized bar strength versus the interaction pa- 
rameter Qe (Elmegreen et al. 1991). The layout is as in Fig. 1151 



and Elmegreen et al. II1991II . however, might be attributed to 
the method for measuring the bar strength. For our purely 
stellar simulations we find a correlation between O and the 
induced bar strength rather than a limiting value of Q. Still, 
the exact correlation between Q and the bar strength - as 
described in Sec. 01 - is likely to be model dependent. 

7.2 Radial redistribution of angular momentum 

Athanassoula l|2()()8|l argued that the redistribution of angu- 
lar momentum is the driving agent for the evolution of an 
isolated barred galaxy. Both her analytical work and her nu- 
merical simulations show that galaxies that have exchanged 
more angular momentum should have a stronger bar, with 
a faster decreasing pattern speed. In our purely stellar mod- 
els, we also find that, whenever bar regeneration occurs, it 
is always accompanied by a considerable loss of angular mo- 
mentum from the disc (sec Fig. 1181 . Indeed, as in the case of 
isolated galaxies (Athariassoula 2003J, there is a tight cor- 
relation between the strength of the bar and the angular 
momentum lost by the inner disc. This argues that angu- 
lar momentum exchange is tightly linked to the bar forma- 
tion process, independent of whether that is spontaneous, 
or driven. 

The actual exchange process, however, is not always 
the same. In isolated discs angular momentum is emitted by 
particles in resonance with the bar in the inner part of the 
disc (within corotation), and, to a lesser extent and if there 
is a considerable bar growth, by non-resonant particles in 
that region. This is absorbed by particles in the outer disc 
and halo. In this case there is only one pattern speed, that 
of the bar, since the spirals, which could in principle have 
a different pattern speed (Tagger et al. 1987; Sygnet et al 
1988), have died away in the early parts of the simulations. 
Thus resonances are well defined. This is certainly not the 
case here. Besides the bar pattern speed, there is the driving 
frequency of the interaction which changes with time. Fur- 
thermore, the forcing from the companion is of comparable 
strength to that of the bar, and the companion itself partici- 
pates actively in the angular momentum exchange. Thus, in 
contrast to the isolated cases, we do not find the bar coro- 



tation radius to separate disc angular momentum emitters 
from disc absorbers. 

The bar regeneration process can be be understood in 
terms of basically the same angular momentum considera- 
tions as described by Athanassoula (12008^ for isolated disc 
galaxies. The bar in an isolated disc can get stronger and 
lose angular momentum basically by four different effects, 
which should be linked to each other. First, particles which 
were on quasi-circular orbits outside the bar get trapped into 
elongated orbits in its outer part and thus the bar becomes 
longer. Secondly, the orbits trapped in the bar could get 
thinner and make the bar thinner, too. Third, more mass 
could be get trapped on periodic orbits in the bar. And 
finally, the bar can of course lose angular momentum by 
slowing down. These effects should also be present in our 
case as well, where the loss of disc angular momentum is 
predominantly driven by the tidal perturbation of the com- 
panion. We will discuss here how much the different effects 
contribute to the regeneration of the stellar bars in our simu- 
lations. Since analytical calculations cannot make any state- 
ments about this, we will base the discussion on the dynam- 
ical properties and orbital structure of our specific models. 

In Fig. 1231 (right-hand column, upper panel) we show 
a direct comparison of the characteristic diagrams of the 
purely stellar models II and II C4p at the end of each run, 
respectively. The characteristic diagram of the isolated and 
the perturbed case look very similar in terms of the layout 
of the main periodic orbits. The main difference, however, 
is that both the xi and the X2 orbits extend to higher en- 
ergies after the interaction. This is especially obvious for 
the part of the xi characteristic around Ej — —0.75, where 
the value of the y intercept increases strongly with Ej. This 
clearly moved towards higher energies. In Fig. 1241 (right- 
hand column) we show some examples of the periodic orbits 
with different Jacobian energies Ej of both the xi and the 
X2 family. We find that the xi orbits become longer and, 
to some smaller extent though, also thinner after the in- 
teraction. This becomes clearer in Fig. 1231 (middle panel), in 
which we plot the axial ratio b/a of the orbits, where a and b 
denote the major and minor axis, respectively, as a function 
of _Ej. The main effect, however, is the lengthening of both 
the xi and the X2 orbits towards higher energies and, tak- 
ing into account the corresponding characteristic diagram, 
also in radial extent. For the xi orbits we also notice some 
lengthening towards lower energies, as well. Owing to these 
effects the bar gets both more centrally concentrated and 
more extended to large radii. 

The thinning of the periodic orbits, which is less pro- 
nounced in our models than the lengthening, can be un- 
derstood by the analytic calculations by Lynden-Bell (1979; 
hereafter LB79) in the context of gradual bar growth in iso- 
lated disc galaxies. As described in LB79, the mean circular 
frequency flc of most disc stars in the central and/or inner 
disc region is much higher than the pattern speed Hp of the 
weak periodic perturbation, which is considered to be bar- 
like in this case. In the frame of reference corotating with 
the perturbation, the fast motion fie — f2p of the stars on 
near-circular orbits is not considerably affected by the weak 
perturbation and therefore the fast action variable Jt (see 
LB79) is approximately constant. As shown by Lynden-Bell, 
angular momentum is removed from the orbits, whose elon- 
gation leads the perturbing bar-like potential. Near-resonant 
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Figure 23. Comparison of models I0C4p (with gas; left-hand side) and II C4p (without gas; right-hand side). The upper panels show 
the characteristic diagram at the end of the run of the interaction model (main panel) and of the corresponding isolated model (subpanel). 
The layout of the diagrams is the same as for, e.g., Fig. |3] Instability regions in the characteristic curves are indicated by dotted lines. In 
the middle panels we show the axial ratios b/a of the xi and X2 orbit families, where a and b are the major and minor axis, respectively. 
The full and dotted lines show the results of the interaction and the isolated model, respectively. The horizontal dashed lines indicate 
axes ratios of 1:2, 1:3, 1:4 and 1:5. The bottom panels show a histogram of the number of stellar disc particles within corotation per 
energy interval, normalized by total number of stellar disc particles within corotation in the corresponding isolated model. The vertical 
dotted lines mark the region in which the X2 orbits is present in the perturbed models. 



orbits in the inner disc region, which reside in the abnormal 
region (see Fig. 2 in LB79), i.e. the region where fii = fie — 1«; 
of the orbit decreases as its angular momentum decreases 
with constant Jf, will align with the perturbation and be- 
come more eccentric, supporting the growth of the bar struc- 
ture. As pointed out by Lynden-Bell, this way the bar shape 
also becomes more eccentric, but the length of the bar will 
not change significantly. The main difference to the scenario 
described in LB79 is that in our case, the perturbation is 
not periodic, since it is due to both the companion and the 
bar, and the angular frequency of the former is a function 
of time. 



We finally check which orbits are important for making 
the bar. In the lower panel of Fig. 1231 we therefore plot the 
number of particles as a function of Ej, normalized by the 
total disc particle number within corotation of the isolated 
model II. This way we get information about how the prin- 
cipal orbits at different energies are populated in our model. 
As can be seen from this plot, the total number of stellar 
particles within the corotation radius increased by less than 
10 per cent. We conclude that the mass of the bar does 
not increase significantly. The main effect we find is that 
particles with energies close to the vertical a;i-branch move 
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IOC4p 




model II (without gas) 




1C4p 




Figure 24. Concentric orbits of the xi (full lines) and X2 (dot- 
ted lines) families. The left- and the right-hand column show the 
models with and without gas, respectively. In the upper panels 
we show the isolated models 10 and II and in the lower panels 
models I0C4p and IlC4p. 



towards higher energies or higher radial extent, taking into 
account the characteristic diagram. 

By tracking individual particle orbits in the correspond- 
ing potential and plotting their average y-intercept versus 
the Jacobian energy Ej in the characteristic diagram, we 
find that the X2 orbits in the purely stellar models are only 
populated by a very small fraction of particles. This result is 
consistent with Fig. 1251 in which we show the surfaces of sec- 
tion of model II C4p at the end of the run for different Ej. 
As can be clearly seen, the xi family dominates phase space 
over a wide range in energy and the X2 orbits are almost not 
present. 



T.3 Influence of the gas 

Interactions which are strong enough to induce a bar in our 
purely stellar models are not sufficiently strong to induce a 
bar in the A'^-body/SPH models. In fact it has not been pos- 
sible to induce a bar in any of our dissipative models, and we 
will here discuss the possible reasons for this. The stellar disc 
component in these simulations loses angular momentum, as 
was found also in the purely A^'-body models (see Fig. Inland 
I18II . The gas in the inner parts, however, gains angular mo- 
mentum, so that the inner parts of the disc may gain a small 
amount of angular momentum (see Fig. |UJ, or if they lose, 
it is considerably less than the corresponding quantity for 
the purely stellar case. Furthermore, as was shown in Fig. [7| 
the interaction is accompanied by a significant inflow of gas 
towards the central disc. These two effects, coupled, pre- 
vent the regeneration of the bar. In this section we discuss 
how the increase of the central mass concentration affects 
the orbital structure of the disc and thus helps prevent the 
regeneration of the bar in these models. 

In Fig. 1231 (left-hand column) we compare models with 
gas, before and after the interaction. In particular, for mod- 
els 10 and I0C4p, we compare the characteristic diagrams 



(upper panel), the axial ratio of the xi and X2 orbits (mid- 
dle panel) and the particle number per energy range (bottom 
panel). In contrast to the models without gas (right column) 
the characteristic diagram of the models with gas extends 
more to lower energies, as would be expected from Fig.|7|and 
Sec. 01 that show clearly that the density distribution is cen- 
trally more concentrated in models with gas. Furthermore, 
the extent of the X2 orbits is considerably larger, and the 
xi orbits have a sizeable instability strip between roughly 
£j = -1.34 and -1.66. 

Information on the axial ratio of the xi and X2 orbits 
and on their extent is given in the middle panels of Fig. 1231 
and in Fig. 1241 We see clearly that the differences between 
the cases before and after the interaction are much smaller 
than for the purely stellar case. The xi orbits are a bit thin- 
ner after the interaction and the X2 orbits have slightly larger 
radial extent, but the differences are small. From Fig. l23l and 
1241 we see that the xi orbits are less elongated than the cor- 
responding orbits in the models without gas. This is partic- 
ularly true at the highest and lowest energies. The extent of 
the elongated orbits is also considerably shorter. Also the X2 
orbits in the central region, i.e. at lower energies, are much 
rounder in the case with gas. 

The number of stellar particles within corotation in 
models with gas (Fig. 1231 bottom panel) have increased by 
roughly 12 per cent compared to the corresponding isolated 
model. Particles move from the intermediate energy region 
to the high energy one, and also, though to a lesser amount, 
to the low energy central region, similar to what is found in 
the models without gas. 

Comparing the SOSs of models with and without 
gas (I0C4p and IlC4p, respectively) after the interaction 
(Fig. 1251 . we note that the area corresponding to xi orbits is 
much smaller in models with gas. On the other hand the area 
corresponding to chaotic motion and the area corresponding 
to X2 orbits is considerably larger. The chaos in the cases 
with gas is due to the instabilities of the xi orbits, discussed 
above. The larger area covered by the X2 orbits is in agree- 
ment with the fact that the X2 characteristic is much more 
extended. To pursue this further we followed in the frozen 
potential the orbits of particles with initial conditions from 
the simulation, as we had already done for the simulations 
without gas. We find a notable difference. Namely there is 
now an indication that, contrary to the purely stellar case, a 
non-negligible fraction of the orbits is trapped around X2 or- 
bits. This difference can be understood as a result of the in- 
duced gas- inflow (see Fig.[7|l, in accordance with the results 
found in numerical simulations of isolated gas-rich barred 
galaxies (e.g., Friedli & Benz 1993; Berentzen et al. 1998). 
Owing to the growing central mass concentration the X2 or- 
bits cover a larger phase-space volume at the expense of the 
xi orbits. 

Thus the fact that bars cannot be regenerated in simula- 
tions with gas can be understood with the help of the many 
differences between the two models, described above. There 
is considerably less angular momentum loss from the inner 
disc material, if this is not a gain. The interaction brings 
considerably less change of the shape and extent of the xi 
orbits, but renders them unstable over a considerable energy 
interval, thus introducing a considerable amount of chaos. 
Finally, due to the increased central concentration, the im- 
portance of the X2 population is considerably increased. The 




Figure 25. Surfaces of section for models I0C4p (witli gas; left-hand column) and IlC4p (without gas; right-hand column) at the end 
of the runs. The values of the Jacobian energy Ej, given in the bottom left corner of each panel, have been chosen such as to give roughly 
the same range in y for each row. In the purely stellar model the phase space is dominated by the xi and 14 orbits at low Jacobian 
energies, while in the models with gas the phase spa<;e is mainly dominated by the X2 and irregular orbits. 
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coupling of these three anti-bar effects prevents the regen- 
eration of a bar component. 



7.4 Properties of regenerated bars 

As a next step we discuss the properties of the regenerated 
bars compared to the ones formed by the bar instabihty in 
isolated discs. Contrary to Miwa & Noguchi (159$) , we find 
no clear evidence for qualitative differences between the two 
types of bars, while we do find clear evidence for similari- 
ties. The orbital differences between the two types of bars, 
discussed in the two previous subsections are quantitative, 
rather than qualitative, since they pertain only to the extent 
of the families and the shape of the corresponding orbits. 
The role of the resonances seems the same. 

Athanassoula II2UU3II found a correlation between the 
bar strength and the pattern speed of bars in isolated disc 
galaxies. In order to compare the properties of the bars 
formed in our models and the ones in isolated discs, we 
first calibrate our model units appropriately, as proposed in 
Athanassoula & Misiriotis ( 2002 I , and then compared them 
to the ones in Athanassoula (120031) . The result is shown in 
Fig. 1261 which shows that the area covered by the driven 
bars is not separated by that covered by the isolated bars. 
This argues strongly for the similarity between the dynami- 
cal properties of the two types of bars, and the difficulty to 
distinguish between them. 

Elmegreen & Elmegreen 1198511 found two difi'erent 
types of major axis surface density profiles, the exponen- 
tial and the flat ones, and Noguchi II199()II argued that this 
distinction could be due to the fact that the first type is 
found in spontaneous bars and the second in driven ones. 
We believe that the difference is due to the bar strength, 
rather than its origin, since Athanassoula & Misiriotis (I2002II 
showed that both types of profiles could be found in isolated 
galaxies: the fiat profiles occurring in the stronger bars and 
the exponential ones in less strong ones. 

Apart from these differences, we find agreement with 
the results of Miwa & Noguchi H1998| l. Thus, in our purely 
stellar simulations we find a roughly linear correlation be- 
tween the pattern speed of the regenerated bar and its 
strength, or, accordingly, the interaction parameter O (see 
Fig. 1171 and 1151 respectively). In other words, stronger bars 
- regenerated by the interaction - have lower pattern speeds. 
This result is in agreement with the simulations of Miwa & 
Noguchi I1998II . who find a similar correlation - less strict 
though - between the pattern speed of the induced bars 
and the companions mass. In both cases, i.e. the formation 
: Noguchi 1998 1 and regeneration of stellar bars (see 
by tidal interactions, the bars always have lower 



l |Miwa . 

Finn 

pattern speeds than the bars formed spontaneously in the 
corresponding isolated models. 

In agreement with Noguchi Ill987t and Gerin et al. 
1199011 we find that tidal bars are not transient, but long- 
lived, as are spontaneous bars. Thus bars which are observed 
at high redshifts ISheth et al. 200311 may have formed by 
galaxy interactions and thus played an important role in 
the evolution of disc galaxies (ISheth et al. 200'3l . 



Figure 26. Bar strength Sb versus pattern speed Qp in units 
used by Athanassoula & Misiriotis (2002). The dots and bullets 
show the results of the isolated models of Athanassoula (2003) 
and our purely stellar interaction models, respectively. 



7.5 Numerical considerations 

The simulations presented in this work have been performed 
with a rather low particle number of the host galaxy, because 
of the relatively large sample of simulations. On the one 
hand the induced numerical noise in general supports the 
bar formation process, while on the other it will heat the 
disc, and thus make the bar regeneration more difficult. We 
trust, however, that a higher number of particles would make 
only quantitative, and not qualitative differences. Thus the 
physical results should remain unchanged. 

The companion galaxy in our simulations has always 
been approximated by a softened point mass. For simula- 
tions in which the companion passes outside the halo of the 
host galaxy, this approximation is well suited. For simula- 
tions in which the companion crosses the halo, however, dy- 
namical friction can be sufficiently strong to change the orbit 
of the companion and finally lead to a merger between the 
two galaxies. 



8 SUMMARY 

In this paper we used numerical simulations to investigate 
the regeneration of a stellar bar by tidal encounters. The 
host galaxy has been chosen to be initially bar unstable and 
forms a large-scale bar during its early evolution. Before the 
interaction with the companion galaxy the bar significantly 
weakened, owing to the bar-driven gas-infiow towards the 
central disc region. For the simulations presented in this 
work we used two different types of host galaxies, i.e. one 
gas-rich disc model and one without gas. The encounters in 
our simulations have been chosen to be prograde with re- 
spect to the rotation in the host disc and co-planar with the 
discs equatorial plane. The mass and the orbit of the com- 
panion have been varied in order to cover a large parameter 
space. 

We found that interactions, which are sufficiently strong 
to regenerate the bar in the purely stellar models, do not lead 
to a regeneration in the dissipative models. The regenerated 
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bars in our simulations are long-living phenomena and by no 
means only transient, i.e. that (regenerated) bars formed by 
interactions may indeed contribute to the number of barred 
spirals at high redshifts and in the local Universe. 

We have shown that the strength of the regenerated 
bars increases with the interaction strength. Owing to the 
tidal perturbation, angular momentum is removed from the 
disc. In fact, the whole disc within its initial cut-off radius 
loses angular momentum, in contrast to what is found for 
isolated discs, where the corotation radius of the bar sep- 
arates disc regions losing and gaining angular momentum. 
The amount of angular momentum removed from the disc 
shows a clear correlation with the interaction strength. Wo 
argued that the main effect of this angular momentum loss 
is a significant extension of the region where bar-supporting 
orbits exist, resulting in a lengthening of the bar. As a fur- 
ther, though somewhat less important, effect, we found also 
a thinning of both the bar-supporting orbits and the bar 
itself. 

The regenerated bars have lower pattern speeds than 
the bars in the corresponding isolated models. Furthermore, 
we found strong correlation between the strength and the 
pattern speed of the bar. This correlation is in very good 
agreement with the correlation found for bars in isolated 
discs. This is one of the pieces of evidence that the regener- 
ated bars are qualitatively similar to those formed in isolated 
discs and thus cannot easily (if at all) bo distinguished from 
them by their dynamical properties. 

In contrast to the purely stellar simulations presented 
in this work, it has not been possible to regenerate the bar in 
our models including gas. We argued that this is due to the 
fact that less angular momentum is lost from the inner disc 
and to the fact that, owing to the interaction, additional gas 
is driven towards the centre of the disc. We concluded that 
the regeneration of stellar bars by galaxy interactions seems 
to be a reasonable mechanism for galaxies containing not 
too much gas, provided the external forcing is sufficiently 
strong. 
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